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Abstract 

Ionization energies for the ground state of berylliumlike ions with nuclear charge numbers in the range 
16 ^ Z ^ 96 are rigorously evaluated. The calculations merge the ab initio QED treatment in the first 
and second orders of the perturbation theory in the fine-structure constant a with the third- and higher- 
order electron-correlation contributions evaluated within the Breit approximation. The nuclear recoil and 
nuclear polarization effects are taken into account. The accuracy of the ionization energies obtained has 
been significantly improved in comparison with previous calculations. 
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I. INTRODUCTION 


During the past four decades, starting from the pioneering works on calculations of one-electron 
self-energy [1] and vacuum polarization [2, 3] corrections, considerable progress has been achieved in 
the theoretical understanding of the calculations of quantum electrodynamics (QED) corrections to 
all orders in aZ, where a is the fine-structure constant and Z is the nuclear charge number. To date, 
state-of-the-art QED calculations include all corrections up to the second order in a, for review 
see Refs. [4-8] and references therein. High-precision experiments to measure the binding and 
transition energies in highly charged ions [9-22] are sensitive to the second-order QED corrections. 
They confirm theoretical predictions to a high accuracy, and thereby provide a unique opportunity 
for tests of QED in the strong Coulomb held. 

Among rigorous QED calculations of the energy levels of highly charged ions, see [23-35] and 
references therein, a particular place is occupied with the works where the ionization energies are 
investigated. The ionization energies allow one to relate masses of ions with different number 
of electrons, and hence are of great interest for the tasks of mass spectrometry [36, 37]. Up to 
date, the high-precision ab initio QED calculations of the ionization energies have been performed 
for hydrogen [32, 38], helium [33, 34], and lithium [35] isoelectronic sequences. There are many 
nonrclativistic and relativistic calculations of the ground-state energies of berylliumlike ions in 
the literature [39-49]. Some of these works somehow include the radiative and nuclear recoil 
corrections. As a rule, only the first-order QED effects are included, the many-electron QED 
effects are treated with the use of one-electron approximations or semiempirical methods. The 
necessity of a more rigorous QED treatment to improve the theoretical accuracy has been pointed 
out in Ref. [49]. 

Recently, in Ref. [50], we have calculated the ground-state binding energies for the even-Z 
berylliumlike ions in the range 18 ^ Z ^ 96. We have also evaluated the ionization energies 
for these ions by subtracting the binding energies of corresponding Li-like ions taken from the 
literature. The main goal of the present investigation is to calculate the ground-state ionization 
potentials for these ions directly and self-consistently including the QED corrections up to the 
second order in a and the electron-correlation effects to all orders in 1/Z. We also extend these 
calculations to all other ions along the beryllium isoelectronic sequence with nuclear charge 16 ^ 
Z ^ 96. 

The paper is organized as follows. The procedure of calculation of the ionization energies is 
described in Sec. II. In Sec. Ill we present our numerical results and compare them with the 
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previous computations. Section IV is reserved for a brief summary. 

The relativistic units (H — c — 1) and the Heaviside charge unit ( a = e 2 /47r,e < 0) are used 
throughout the paper. 

II. BASIC FORMULAS 

The Furry picture is the standard starting point for investigation of the properties of heavy 
few-electron ions within QED. In this picture interaction of electrons with the Coulomb field of 
the nucleus V^uci is taken into account to all orders in aZ from the very beginning: 

[-icx. • V + f3m + Kuci(r)] i’n(r) = e n ^ n {r). (1) 

The interelectronic interaction and QED effects have to be accounted for by perturbation theory. 
A convenient approach to construct the QED perturbation series is the two-time Green function 
(TTGF) method [51]. 

Naturally, such an approach works well when the interelectronic interaction is small in compar¬ 
ison with the binding energies of all electrons of the system. Otherwise, the convergence of the 
perturbation series is often slow. It is possible to accelerate the convergence of the perturbation 
theory by transition to the extended Furry picture. This implies the replacement of the nuclear 
potential in Eq. (1) with an effective potential: 

Kuci(r) -t V eS (r) = U nU ci(r) + V SCT (r). (2) 

The local screening potential 14 cr (r) in Eq. (2) models the screening of the nuclear potential by 
the other electrons partly accounting for the interelectronic interaction in the zeroth-order Hamil¬ 
tonian. The interaction with the potential 5V(r) = — V scr (r) is to be accounted for perturbatively 
to avoid the double counting of the screening effects. Furthermore, the calculation of the ls 2 2s 2 
state energy starting from the pure Coulomb field as the zeroth order approximation is difficult, 
because of the quasidegeneracy with the ls 2 (2p!/ 2 ) 2 state. The application of some screening po¬ 
tentials can remove this degeneracy and simplify calculations significantly. The extended Furry 
picture was successfully employed in QED calculations of the energy levels [28, 30, 31, 35, 52-55], 
the hyperfine splitting [56-61], and the g -factor [62-64], 

In this paper we apply four different types of the screening potential. All these potentials are 
spherically symmetric by construction. Our first and simplest choice is the core-Hartree (CH) 
potential induced by the Is 2 closed shell. It can be obtained from the radial charge density of two 
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Is electrons: 


r°° i 

Vch (r) = a dr' — pcn(r% (3) 

Jo r > 

/»oo 

Pch(t) = 2 [G? a (r) + F 2 s (r)] , / pca(r)dr = 2, (4) 

Jo 

here G/r and F/r are large and small radial components of the Dirac wave function. 

Three other types of the screening potential were constructed for both three-electron (ls 2 2s) 
and four-electron (ls 2 2s 2 ) configurations, so that the total number of the screening potentials 
employed is equal to seven. In what follows we label these potentials with the indices 3 and 4, 
depending on the electronic configuration used. The first type of the screening potential generated 
for both configurations is the local Dirac-Fock (LDF) potential [65]. This potential is constructed 
by inversion of the radial Dirac equation with the radial wave functions obtained in the Dirac-Fock 
approximation. Two other screening potentials arise from the density-functional theory. In terms 
of the total radial charge density pt of all electrons in the configuration the Kohn-Sham (KS) 
potential reads as follows: 

Vjss( r ) = a f o dr '^rPt( r ') - “ (3^ rpt ( r )) ’ ( 5 ) 

/♦oo 

Pt(r) = p C n(r) + (N - 2) [G 2 s (r) + F 2 2 (r)] , / p t (r)dr = N, (6) 

Jo 

where N = 3 or 4 depending on the electron configuration. We introduce the Latter correction 
[66] to improve the asymptotic behavior of the KS potentials at large r. Finally, the last type of 
the screening potential generated for ls 2 2s and ls 2 2s 2 configurations is the Perdew-Zunger (PZ) 
potential [67]. Summarizing the description of the screening potentials used in the present work 
it is worth noting that they have different asymptotic behavior. The potentials Vch, Pks 3 , Vldf 3 
and Vpz 3 behave like 2 a/r at large distances r, whereas the potentials Vks 4 , Vldf 4 , Vpz 4 behave 
like 3 a/r. In other words, if we set N SCI = N — 1 for the LDF, KS and PZ potentials, and N scr = 2 
for the CH potential, then all screening potentials behave like N scr a/r. In this case eN scr has the 
meaning of the total charge of the screening cloud. 

The calculation of the ionization potentials of berylliumlike ions in the present work follows 
in general the scheme which was used in Ref. [50] for the evaluation of the binding energies. 
Nevertheless, there are some modifications in the calculation procedure. Below, we briefly discuss 
the main stages of the calculations. 

The evaluation of the ionization potentials can be conveniently divided into several steps. At 
first stage it is necessary to solve Eq. (1) for the state under consideration with an effective 
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potential (2). Moreover, in the bound-state QED calculations one needs to have a quasicomplete 
basis set of the Dirac equation solutions for representation of the Dirac-Coulomb Green function 


G(u,r l ,r 2 ) = Y 

n 


^n(n)^(r 2 ) 

OJ - £ n (l - iO)' 


(7) 


The set of one-electron wave functions was evaluated employing the dual kinetic balance (DKB) 
method [68] with the basis functions constructed from the B-splines [69]. 

The application of the screening potentials allows one to partly take into account the interaction 
between electrons. We have to consider the remaining interelectronic interaction by perturbation 
theory. In Fig. 1 the corresponding set of first- and second-order Feynman diagrams is depicted. 
The double line represents the electron propagator in the effective potential. The circle with a cross 
corresponds to the screening potential counterterm. The ionization energy of the 2s electron can be 
obtained by subtracting the binding energy of the Li-like ion from the binding energy of the Be-like 
ion. Diagrams containing only the interaction between the Is electrons contribute to both binding 
energies and cancel each other in the difference. For this reason, here and further throughout the 
paper we deal only with those diagrams, where one or two of the incoming (outgoing) electron 
lines correspond to the 2s electrons. The calculation formulas for the diagrams (a)-(d) in Fig. 1 
can be found, e.g., in Refs. [26, 70]. As it was noted in Ref. [50], only a minor modification of the 
standard procedures is necessary to perform the calculation of the 2s 2 interaction. The formulas 
for the corrections corresponding to the diagrams (e)-(g) can be easily obtained with the use of 
the TTGF method: 
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( 8 ) 

(9) 


( 10 ) 


where a = 2 s, b = Is, V ab = {a\5V\b) = (a| -V sa \b), I abc d(u) = (ab\I (u)\cd), /(ce) = e 2 a l ‘a l 'D ltt ,(i t j), 
D^ w (cu) is the photon propagator, I' abcd (u >) = (ab| ■£l(uj)\cd),I ab . cd = (ab\I(A bd )\cd)-(ba\I(A ad )\cd), 
and A ab = e a - £ b . 
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FIG. 1. The interelectronic-interaction diagrams. The double line denotes the electron propagator in the 
effective potential (2). The wavy line corresponds to the photon propagator. The symbol <8> represents 
the local screening potential counterterm. 


Next, we have to take into account the interelectronic-interaction corrections of the third and 
higher orders. We have evaluated these contributions within the Breit approximation. At first, we 
have calculated the binding energies of the ls 2 2s and ls 2 2s 2 states solving the Dirac-Coulomb-Breit 
equation by means of the configuration-interaction Dirac-Fock-Sturnr method (CI-DFS) [71, 72], 
Then, the contributions of third and higher orders to these binding energies, E{' 2 ’ 2 and 2 , 
were extracted from the total CI-DFS energies using the numerical procedure described, e.g., in 
Refs. [30, 31]. Finally, we took the difference E^ 2s2 — E^^ to obtain the desired contribution 
E\nt Breit to the ionization potential of the 2s electron. To estimate the accuracy of our numerical 
procedure we have compared the first- and second-order interelectronic-interaction contributions to 
the ionization potential, E^ Breit and ¥nt,Breit’ calculated within CI-DFS method with analogous 
quantities evaluated with the use of our code for the QED calculations. For this reason the diagrams 
(a), (b), (d)-(g) in Fig. 1 have been evaluated in the Coulomb gauge neglecting the contribution of 
the negative-energy continuum and energy dependence of the interelectronic-interaction operator. 
We have found a very good agreement between those two independent approaches for all types of 
the screening potentials. 

Next, one should consider the contribution of the radiative corrections. Figure 2 displays all 
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relevant first- and second-order QED diagrams with the exception of one-electron two-loop graphs, 
which we will discuss below. The diagrams from the first line are known as the self-energy (SE) 
diagrams, while the diagrams in the second line are referred to as the vacuum polarization (VP) 
diagrams. The formal expressions for these contributions can be easily obtained within the TTGF 
method, they are presented, e.g., in Ref. [30]. The formal expressions suffer from ultraviolet 
divergences. The way for cancellation of the divergences in the SE diagrams is discussed in detail 
in Refs. [73, 74], in the present work we follow the renormalization scheme described there. 

The VP contributions are conveniently divided into the Uchling and Wichmann-Kroll parts. 
The renormalized expression for the Uchling potential is well known: 


Vijehl(?') 


- aZ ^l d/4xrV ~ l(r ')/ 

{exp(— 2m\r — r'\t ) — exp(—2 m(r + r')t)} 
Arrirt 


(ii) 


where p n uci(^) is the nuclear charge density, normalized to unity. To account for the screening effect 
on the Uchling potential in the diagrams (f), (i)-(j), one should replace Zp nucl (r) by Zp nucl (r) — 
N s C rPscr(r), where p scr (r) is the charge density of the electron screening cloud, normalized to unity. 
In order to calculate the contribution of the counterterm diagram (h) one has to replace Zp nuc ] (r) 
by -VscrPscr (j") hr Eq. (11). The renormalized Uehling operator which corresponds to the diagram 
(g) has the form 


Iuehi(£,r,r') 


a^a^ 2 a f°° 
| r -r' | 3 7r J 1 
x exp (—a/ (2 mt) 2 


dt (* + Jfi) 

- e 2 \r - r'|), 


V^i 

t 2 


( 12 ) 


where £ is the energy of the transmitted photon. 

The Wichmann-Kroll contribution was calculated for the diagrams (f), (i) and (j) using the 
expression: 


2a 


zboo 


IVkM = - £ 
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x / duj I dy y 2 dz z 2 -7-- 

Jo Jo Jo max(r, y) 

2 

xUeff(^) Re iK k ( iuj iy’ z ){ G t(, iuj ’yi z ) 

i,k =1 




(13) 


where F f \ k and Gf' : are the radial components of the partial contributions to the free- and bound- 
electron Green function. It turns out that calculation of Eq. (13) using the spectral decomposition 
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FIG. 2. First- and second-order QED diagrams (excluding the one-electron two-loop diagrams). The 
notations are the same as in Fig. 1. 


(7) with the B-spline finite basis set is rather problematic. Therefore, we have evaluated the 
radial Green functions for all screening potentials numerically, solving a corresponding system of 
differential equations. The Wichmann-Kroll contribution to the diagram (g) is relatively small 
[75, 76] and has been neglected, together with the related contribution to the diagram (h). 

The consideration of the two-loop one-electron contributions completes the rigorous QED treat¬ 
ment in the second order in a. Evaluation of these corrections to all orders in aZ is a very de¬ 
manding problem which has not yet been solved completely. The most remarkable progress in this 
area was achieved by Yerokhin et al. [29, 77]. In Ref. [77] the complete gauge-invariant set of the 
two-loop self-energy diagrams was calculated for high-Z ions. The so-called SEVP, VP VP, and 
S(VP)E sets were tabulated in Ref. [29] using the free-loop approximation if the rigorous calculation 
was not possible. We take the two-loop one-electron corrections from Refs. [29, 38, 77]. 

All corrections considered so far correspond to the approximation of an infinite nuclear mass. 
Nevertheless, high precision calculations of energy levels in highly charged ions must include the 
nuclear recoil corrections. The full relativistic theory of the recoil effect can be formulated only in 
the framework of QED. Such a theory to the first order in m/M (M is the nuclear mass) and to 
all orders in aZ was worked out in Refs. [78, 79] (see also Refs. [80, 81] and references therein). 
Within the Breit approximation the following many-body Hamiltonian describes the contribution 



























of the recoil effect to the binding energy [78, 79, 82] 


hi 

Evaluating the expectation values of the relativistic operator (14) with the CI-DFS wave functions 
of 1s 2 2s and ls 2 2s 2 states and taking the difference we obtain the nuclear recoil contribution to 
the ionization potential of the 2s electron in the Breit approximation. 

The nuclear recoil effects which are beyond the Breit approximation are termed as the QED 
recoil effects. In the present work we calculate these corrections in the zeroth order in 1/Z. To 
this order only one-electron terms contribute. The two-electron terms vanish, since the Is and 2s 
electrons have the same parity [83]. For the point-charge nucleus, the QED recoil correction to 
first order in m/M, to zeroth order in a, and to all orders in aZ is given by [78, 79]: 


Pi ' Pj 


aZ 


OLj + 


ctj ■ r, r,- 


Pj 


(14) 


A E 
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2i tM J _<*, 

xG(w + £ a ) ( 
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UJ 


Ip, v ] A 

u> + iO / 


D(co) + 


IP, V] \ 

u + iO ) 



(15) 


Here p is the momentum operator, D m (u> ) = —AnaZaiDi m {uj), and Di m is the transverse part 
of the photon propagator in the Coulomb gauge. The nuclear size corrections to the QED recoil 
effect can be partly accounted for by replacement of the potential V, the Dirac energy e a , the 
wave function |a), and the Green function G in Eq. (15) with the corresponding quantities for the 
extended nucleus [80]. Such calculations were carried out in Refs. [84, 85] for Is and 2s states of 
H-like ions. In the present work we have performed the calculations of the QED nuclear recoil 
corrections for the effective potential (2) in order to take into account the screening effects. 

Finally, one has to consider the nuclear polarization effects. These effects are related to the 
intrinsic nuclear dynamics. We have taken into account this correction for all ions using the results 
of Refs. [86-90] and the corresponding prescriptions from Ref. [38]. 

We note that the numerical procedures are checked by performing the calculations of the QED 
corrections in two gauges, the Feynman and the Coulomb ones. A good agreement between both 
calculations is found. The Fermi model with a thickness parameter of 2.3 fin is used to describe 
the nuclear charge distribution. The nuclear radii are taken from Refs. [91]. In the cases of Z with 
no experimental radii available, we use the root-mean-square (rms) radii obtained by employing 
the approximate formula from Ref. [32] and prescribe a 1% uncertainty to them. 
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TABLE I. Individual contributions to the ionization potential of the 2s electron in berylliumlike calcium 
(in eV). See text for details. 


Contribution 

CH 

LDF3 

KS3 

PZ3 

LDF4 

KS4 

PZ4 

#(°) 

■^Dirac 

-1147.45480 

-1149.96374 

-1154.32591 

-1157.78287 

-1074.07837 

-1066.78185 

-1079.51335 

771 ( 1 ) 

Ant 

65.89672 

68.74389 

73.56661 

77.27720 

-7.77584 

-14.87760 

-2.02387 

£l( 2 ) 

^int, Breit 

-9.77636 

-9.37458 

-9.52861 

-12.17568 

-7.48702 

-7.08279 

-9.28264 

Unt,QED 

0.00057 

0.00055 

0.00053 

0.00054 

0.00055 

0.00051 

0.00054 

^O 3 ) 

^int, Breit 

3.88459 

3.14463 

2.83836 

5.23177 

1.89125 

1.29221 

3.37002 

Ant, total 

-1087.44927 

-1087.44925 

-1087.44902 

-1087.44903 

-1087.44944 

-1087.44952 

-1087.44931 

e hi 

^SE 

0.19895 

0.19917 

0.20072 

0.20187 

0.19207 

0.18907 

0.19398 

Z7l(l) 

-0.01388 

-0.01388 

-0.01399 

-0.01408 

-0.01340 

-0.01320 

-0.01354 

^ScrSE 

-0.02115 

-0.02142 

-0.02314 

-0.02439 

-0.01397 

-0.01073 

-0.01600 

^ScrVP 

0.00146 

0.00146 

0.00158 

0.00168 

0.00096 

0.00074 

0.00111 

Z7l(2) 

- C/ 21oop 

-0.00016 

-0.00016 

-0.00016 

-0.00016 

-0.00016 

-0.00016 

-0.00016 

-^QED,total 

0.16521 

0.16516 

0.16501 

0.16491 

0.16549 

0.16573 

0.16538 

-Aec,Breit 

0.01415 

0.01415 

0.01415 

0.01415 

0.01415 

0.01415 

0.01415 

-Aec,QED 

0.00006 

0.00006 

0.00006 

0.00006 

0.00006 

0.00006 

0.00006 

-Aotal 

-1087.26985 

-1087.26988 

-1087.26980 

-1087.26990 

-1087.26974 

-1087.26958 

-1087.26972 


III. NUMERICAL RESULTS AND DISCUSSIONS 


The individual contributions to the 2s electron ionization potentials of berylliumlike calcium, 
xenon, and uranium evaluated for various screening potentials are presented in Tables I—III, re¬ 
spectively. The first line contains the ionization potentials obtained by solving Eq. (1) for the 
corresponding effective potential with the finite nuclear size effect included, Ej-^. ac = e 2s — me 2 . 
The nuclear deformation correction was added for uranium ion in accordance with the results of 
Ref. [92], E^J stands for the contribution of the first-order diagrams in Fig. 1 [diagrams (a) and 
(e)]. The diagram (a) was evaluated keeping the energy dependence of the photon propagator in 
the exchange part, i.e., in the framework of QED. The contribution of the second-order diagrams 
evaluated in the Breit approximation is given in the third line. As it was mentioned in Sec. II, 
here the Breit approximation implies the calculations performed in the Coulomb gauge at zero 
energy transfer and neglecting the contribution of the negative-energy continuum. We note that 
this approach to the Breit approximation differs, e.g., from the one used by Yerokhin et al. in 
Ref. [55], where the exchange by two Breit photons was considered as belonging to the higher order 
corrections. Furthermore, the influence of the negative energy continuum was partly accounted 
for in Ref. [55] in the contribution considered. In the fourth row we give the QED correction 
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TABLE II. Individual contributions to the ionization potential of the 2s electron in berylliumlike xenon 
(in eV). See text for details. 


Contribution 

CH 

LDF3 

KS3 

PZ3 

LDF4 

KS4 

PZ4 

E 

■^Dirac 

-9786.2460 

-9796.3872 

-9809.3421 

-9821.1550 

-9569.5159 

-9548.0653 

-9587.1179 

E (1) t 

198.9440 

209.5581 

223.0184 

235.1590 

-18.0970 

-39.3573 

-0.1378 

_g( 2 ) 

^int, Breit 

-13.6879 

-13.2832 

-13.4699 

-16.6451 

-11.1410 

-10.7201 

-13.2829 

z?( 2 ) 

■ C/ int,QED 

0.0310 

0.0302 

0.0292 

0.0298 

0.0306 

0.0299 

0.0305 

E ’O 3 ) 

^int, Breit 

3.9484 

3.0718 

2.7540 

5.6013 

1.7125 

1.1030 

3.4975 

-^int, total 

-9597.0104 

-9597.0103 

-9597.0105 

-9597.0101 

-9597.0108 

-9597.0098 

-9597.0106 

-^SE 

7.2387 

7.2471 

7.2704 

7.2871 

7.1549 

7.1218 

7.1828 

771 ( 1 ) 

-^VP 

-0.9372 

-0.9382 

-0.9416 

-0.9438 

-0.9265 

-0.9229 

-0.9305 

E ( 2 ) 

-^ScrSE 

-0.2890 

-0.2979 

-0.3221 

-0.3395 

-0.2039 

-0.1700 

-0.2327 

E ( 2 ) 

^ScrVP 

0.0380 

0.0391 

0.0427 

0.0450 

0.0272 

0.0235 

0.0313 

E^ 2 ) 

■^21oop 

-0.0165 

-0.0165 

-0.0165 

-0.0165 

-0.0165 

-0.0165 

-0.0165 

-^QED, total 

6.0340 

6.0336 

6.0328 

6.0322 

6.0351 

6.0358 

6.0343 

-^rec, Breit 

0.0386 

0.0386 

0.0386 

0.0386 

0.0386 

0.0386 

0.0386 

-®rec,QED 

0.0031 

0.0031 

0.0032 

0.0032 

0.0031 

0.0031 

0.0031 

-^'nucl.pol. 

- 0.0002 

- 0.0002 

- 0.0002 

- 0.0002 

- 0.0002 

- 0.0002 

- 0.0002 

-®total 

-9590.9349 

-9590.9352 

-9590.9362 

-9590.9363 

-9590.9342 

-9590.9325 

-9590.9347 


.EjntQED t° tl ie third line, E-J Rr ei t • It is calculated as the difference between the contribution of 
the second-order diagrams in Fig. 1 calculated within the rigorous QED approach and within the 
Breit approximation. The interelectronic-interaction corrections EE t g reit of the third and higher 
orders evaluated with the CI-DFS method are shown in the fifth row. In the sixth line the sum 
Ei n t ,total °f all the previous terms (the rows from first to fifth) is presented. From Tables I III 
it can be seen that for all seven screening potentials the Eint,total values are in a good agreement 
with each other. The next four rows contain the contributions of the QED diagrams from Fig. 2. 
The lines labeled as Egg and Eyp correspond to the first-order SE and VP corrections, whereas 
the lines Eg c ( SE and Eg c ( VP present the related SE and VP two-electron corrections. The two-loop 
one-electron corrections for the 2s electron are collected together in the eleventh line, E^ op . The 
sum of all QED contributions (the rows from seventh to eleventh) is given in the line labeled with 
Eqed. total- Again, it is seen that the results of the calculations are in a good agreement for all 
screening potentials. One can check that this holds independently for the SE and VP corrections. 
The next two rows contain the nuclear recoil contribution evaluated within the Breit approximation 
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TABLE III. Individual contributions to the ionization potential of the 2s electron in berylliumlike uranium 
(in eV). See text for details. 


Contribution 

CH 

LDF3 

KS3 

PZ3 

LDF4 

KS4 

PZ4 

E 

■^Dirac 

-32827.813 

-32858.372 

-32886.232 

-32913.911 

-32407.665 

-32369.166 

-32445.215 

E (1) t 

406.898 

438.333 

466.938 

495.120 

-13.679 

-51.938 

24.367 

_g( 2 ) 

■^intjBreit 

-18.134 

-17.775 

-18.129 

-23.634 

-14.308 

-13.906 

-17.824 

z?( 2 ) 

■ C/ int,QED 

0.275 

0.273 

0.270 

0.272 

0.278 

0.277 

0.277 

E ’O 3 ) 

'Ant,Breit 

6.398 

5.165 

4.778 

9.778 

2.997 

2.352 

6.018 

-^int, total 

-32432.376 

-32432.376 

-32432.375 

-32432.375 

-32432.378 

-32432.382 

-32432.377 

-^SE 

62.677 

62.770 

62.923 

63.026 

62.280 

62.185 

62.455 

771 ( 1 ) 

-^VP 

-14.951 

-14.980 

-15.029 

-15.049 

-14.863 

-14.855 

-14.910 

E ( 2 ) 

-^ScrSE 

-1.746 

-1.843 

- 2.002 

-2.109 

-1.346 

-1.250 

-1.526 

E ( 2 ) 

^ScrVP 

0.423 

0.454 

0.506 

0.526 

0.335 

0.329 

0.383 

E^ 2 ) 

■^21oop 

-0.244 

-0.244 

-0.244 

-0.244 

-0.244 

-0.244 

-0.244 

-^QED, total 

46.158 

46.156 

46.154 

46.150 

46.162 

46.164 

46.158 

-^rec,Breit 

0.066 

0.066 

0.066 

0.066 

0.066 

0.066 

0.066 

-®rec,QED 

0.047 

0.047 

0.047 

0.047 

0.046 

0.046 

0.046 

-^'nucl.pol. 

-0.036 

-0.036 

-0.036 

-0.036 

-0.036 

-0.036 

-0.036 

-®total 

-32386.141 

-32386.142 

-32386.144 

-32386.148 

-32386.139 

-32386.141 

-32386.142 


and the QED recoil correction, respectively. For uranium and xenon we also give the contribution 
of the nuclear polarization effect in the line i? nuc i. po i.. Finally, in the last line we present the total 
results for the 2s electron ionization energies. From Tables I III it can be seen, that the total values 
of the ionization potentials are almost independent of the type of the screening potential used in 
the calculations. For this reason, for all other ions in the range 16 ^ Z ^ 96 the calculations have 
been performed for two screening potentials only, namely for 14 df 3 and Lldf 4 - We have chosen 
these two potentials in order to control the accuracy of the calculations along the isoelectronic 
sequence. 

Table IV displays the ionization potentials for all berylliumlike ions with Z = 16 — 96. For cal¬ 
cium, xenon and uranium we give the averages of the total values evaluated for the seven screening 
potentials. For the other ions the ionization potentials were obtained by averaging the results 
of calculations with the LDF3 and LDF4 potentials. The uncertainties given in the parentheses 
were obtained by summing quadratically the uncertainty due to the nuclear size effect, the un¬ 
certainty of the CI-DFS calculation, and the uncertainties due to uncalculated higher-order QED 
contributions. For uranium ion the nuclear size uncertainty was estimated in accordance with 
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Ref. [92] taking into account the nuclear deformation correction. For the other ions we estimated 
this uncertainty by adding quadratically two contributions. The first one was calculated by vary¬ 
ing the rms nuclear radius within its error bar. The second one rather conservatively estimates 
the uncertainty in the nuclear charge distribution. It was obtained by varying the distribution 
model from the Fermi model to the homogeneously charged sphere model. The uncertainty due 
to uncalculated QED corrections to the interelectronic-interaction contributions can be conserva- 
tively estimated as the product of the term E^ t g reit and the doubled ratio of the second-order 
interelectronic-interaction QED correction to the corresponding contribution calculated within the 
Breit approximation, E int [ q RD /E int [ Breit . The contribution of the higher-order screened QED dia¬ 
grams can be estimated by multiplying the second-order QED term by the factor 2 jZ. On the 
other hand, the total results should be fully independent of the zeroth-order approximation po¬ 
tential, provided the interelectronic-interaction and QED corrections are calculated to all orders. 
Therefore, the discrepancy between the calculations with different potentials may serve as an esti¬ 
mation of the uncalculated QED contributions. We have used both estimations of the uncalculated 
QED corrections and summed them quadratically. Finally, we have estimated the uncertainty due 
to uncalculated two-loop one-electron QED contributions according to Ref. [29]. For \ow-Z ions 
the main source of the uncertainty is the CI-DFS calculation. For liigh-Z ions the nuclear size 
uncertainty and the uncertainties due to uncalculated QED corrections start to play an important 
role. 


In Table IV we compare our results for the ionization potentials of berylliumlike ions with 
theoretical predictions by other authors. It can be seen that, as a rule, our ionization potentials are 
in a good agreement with the results of the previous calculations, but have much higher accuracy. 
This is not surprising since, as was noted in Sec. I, in all previous evaluations of the ionization 
potentials of Be-likc ions the QED effects were included either semiempirically or within some 
one-electron approximations. In our previous work [50] we tabulated the ionization potentials for 
even-Z ions in the range 18 ^ Z ^ 96. The ionization potentials were obtained as the differences 
between the binding energies of Be-like ions calculated directly in Ref. [50] and the binding energies 
of Li-like ions compiled from the previous tabulations. In the present work the ionization potentials 
were calculated directly and self-consistently. This resulted in a three-times improvement of the 
theoretical precision for heavy ions and in a two-times improvement for low- and middle-Z ions 
compared to our previous calculations [50]. 
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TABLE IV. Ionization potentials (in eV) for berylliumlike ions with Z = 16 — 96. 


Nucleus This work 

?§S -652.135(20) 


?SfCl -750.465(20) 


f°Ar -855.750(20) 


fgK -968.010(20) 


2oCa -1087.270(20) 


45 

21 


Sc 


iTi 


51 

23 


V 


52 

24 


Cr 


25 Mn 


-1213.555(20) 


-1346.889(20) 


-1487.301(20) 


-1634.820(20) 


-1789.478(20) 


Other works 
—651.96(12) a 
-652.20 b 
—652.1923 c 
—652.3391 d 
—750.23(19) a 
—750.54 b 
—750.5212 c 
—750.7090 d 
—855.47(27) a 
—855.82 b 
—855.8045 c 
—856.0432 d 
—967.66(35) a 
—968.08 b 
—968.0597 c 

— 1086.85(40) a 

— 1087.3 b 

— 1087.311 c 

— 1087.44 e 

— 1213.06(45) a 
-1213.6 b 
-1213.583 c 

— 1346.33(47) a 

— 1347.0 b 

— 1346.899 c 

— 1486.67(52) a 

— 1487.4 b 
—1487.285 c 

— 1634.11(55) a 
—1634.9 b 

— 1634.769 c 

—1788.70(62) a 
—1789.6 b 
—1789.381 c 


14 



TABLE IV. (Continued.) 


Nucleus 

This work 

Other works 

iFe 

-1951.307(21) 

— 1950.4(1.8) a 

—1951.4 b 

-1951.10 f 

f 7 Co 

-2120.343(22) 

-2119.4(1.9) a 

—2120.5 b 

—2120.14 f 

£ 

0000 

LO(M 

-2296.621(23) 

-2295.6(2.l) a 

—2296.7 b 

-2296.42 f 

—2298.27 e 

iCu 

-2480.182(23) 

-2479.1(2.2) a 

—2480.3 b 

—2479.98 f 

30 ^ 11 

-2671.064(24) 

—2669.9(2.5) a 

-2671.2 b 

—2670.84 f 

—2673.10 6 

3?Ga 

-2869.311(25) 

—2868.1(2.8) a 

—2869.5 b 

-2869.08 f 

32 Ge 

-3074.967(25) 

—3073.6(3.0) a 

—3075.1 b 

-3074.73 f 

slAs 

-3288.077(26) 

—3286.6(3.2) a 

—3288.2 b 

—3287.84 f 

CO 00 

c n 

CD 

-3508.691(27) 

—3507.1(3.5) a 

—3508.9 b 

—3508.45 f 

35 Br 

-3736.857(27) 

—3735.3(3.7) a 


—3737.0 b 
—3736.65 f 
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TABLE IV. (Continued.) 


Nucleus 

This work 

Other works 

36 K- 1 ' 

-3972.631(28) 

—3970.8(4.0) a 

—3972.8 b 

—3972.45 f 

—3975.7297 d 

iRb 

-4216.065(28) 

—4214.2(4.2) a 

—4216.3 b 

IfSr 

-4467.217(29) 

—4465.3(4.5) a 

—4467.4 b 

89 v 

39 1 

-4726.144(30) 

—4724.0(4.7) a 

—4726.4 b 

907 r 

40 LT 

-4992.909(30) 

—4990.7(5.0) a 

—4993.1 b 

41 Nb 

-5267.575(31) 

—5265.2(5.2) a 

—5267.8 b 

42 Mo 

-5550.207(31) 

—5547.8(5.5) a 

—5550.5 b 

97 T 

43 1C 

-5840.876(32) 

—5838.3(5.6) a 

-5841. l b 

44 2 Ru 

-6139.651(33) 

—6136.8(5.8) a 

—6139.9 b 

ifRh 

-6446.608(33) 

—6443.7(6.l) a 

—6446.9 b 

|° 6 Pd 

-6761.822(34) 

—6758.8(6.3) a 

-6762.l b 

47 7 Ag 

-7085.373(34) 

—7082.1(6.6) a 

—7085.7 b 

114 Crl 

48 

-7417.341(34) 

—7413.9(6.8) a 

—7417.7 b 

ii5 Tn 

49 111 

-7757.815(34) 

—7754.2(7.1) a 

—7758.2 b 

50°Sn 

-8106.879(35) 

—8103.1(7.3) a 

—8107.2 b 

^Sb 

-8464.627(35) 

—8455(4) g 

—8465.0 b 
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TABLE IV. (Continued.) 


Nucleus This work Other works 


52° Te 

-8831.151(35) 

-8821(4)8 

—8831.5 b 

Ifl 

-9206.555(36) 

-9196(4)8 

—9207.0 b 

54 2 Xe 

-9590.935(37) 

-9581(4)8 

—9591.4 b 

IfCs 

-9984.401(37) 

-9974(4)8 

—9984.9 b 

—10002.68 e 

ifBa 

-10387.058(38) 

-10376(4)8 

—10388 b 

139 T a 

57 - La 

-10799.024(39) 

-10789(5)8 

—10800 b 

140 

58 

-11220.414(40) 

-11210(5)8 

—11221 b 

141 p_ 

59 rr 

-11651.351(41) 

-11641(5)8 

—11652 b 

142 Nd 

60 iNQ 

-12091.960(42) 

-12082(5)8 

—12092 b 

If Pm 

-12542.365(47) 

-12532(6)8 

152 a 

62 

-13002.682(45) 

-12992(6)8 

153 ■p.. 

63 

-13473.107(46) 

-13462(7)8 

158PU 

64 

-13953.751(48) 

-13943(10)8 

i 9 Tb 

-14444.844(97) 

-14434(12)8 

66 4 Dy 

-14946.339(52) 

-14936(15)8 

67 5 H0 

-15458.633(58) 

-15448(19)8 

68 6 El' 

-15981.752(56) 

-15971(24)8 

gg 9 Tm 

-16515.973(59) 

-16505(30)8 

70 4 Yb 

-17061.335(61) 

-17050(40)8 

7i 5 Lu 

-17618.130(71) 

-17607(40)8 

72° Hf 

-18186.626(67) 

-18176(50)8 

73'Ta 

-18766.923(70) 

-18756(50)8 

^ 4 W 

-19359.249(74) 

-19362.5(3.1) 

-19348(50)8 
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TABLE IV. (Continued.) 


Nucleus 

This work 

Other works 

ifRe 

-19963.856(81) 

—19953(60) g 

7(POs 

-20580.902(81) 

—20570(70) g 

193t,. 

77 11 

-21210.77(22) 

—21200(90) g 

^Pt 

-21853.579(89) 

—21843(90) g 

79 7 Au 

-22509.658(95) 

—22498(100) g 

202 tr„ 

80 P-S 

-23179.23(10) 

—23168(110) g 

IfTl 

-23862.66(10) 

—23852(120) g 

82 8 Pb 

-24560.15(11) 

—24548(120) g 

i 9 Bi 

-25272.10(12) 

—25260(140) g 

il°Po 

-25998.67(14) 

—25988(150) g 

is 5 A t 

-26740.49(28) 

—26729(160) g 

220 r> 

86 nn 

-27497.35(17) 

—27486(170) g 

87 3 F' r 

-28270.33(18) 

—28259(190) g 

88 6 Ra 

-29059.47(24) 

—29048(200) g 

89 7A C 

-29865.71(41) 

—29854(220) g 

232 rpi 

90 ln 

-30687.91(20) 

—30677(240) g 

231 Pa 

91 ra 

-31528.98(50) 

—31516(250) g 

ifu 

-32386.14(20) 

—32374(300) g 

93 7 N P 

-33263.97(61) 

—33252(300) g 

94°PU 

-34158.55(48) 

—34147(300) g 

if Am 

-35073.44(27) 

—35062(300) g 

ifCm 

-36009.57(37) 

—35996(400) g 


a Biemont et al. [40] 
b Gu [43]. 
c Chung et al. [39]. 
d Pathak et al. [48]. 
e Chaudhuri et al. [41]. 
f Yong-Qiang et al. [46] 

g Rodrigues et al. [42] with the uncertainty prescribed by NIST [93]. 
h Kramida and Reader [45]. 
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IV. SUMMARY 


To summarize, we have performed the ab initio QED calculations of the ground-state ionization 
energies for all berylliumlike ions in the range 16 ^ Z ^ 96. Our numerical approach merges 
the rigorous QED evaluations in the first and second orders of the perturbation theory with the 
calculations of the higher-order electron-correlation contributions within the large-scale CI-DFS 
method. As the result, we have obtained the most precise theoretical predictions for the ionization 
potentials in Be-like ions. The achieved accuracy allows to probe the QED corrections in the 
ionization energies of berylliumlike ions and thus provides a prerequisite for tests of QED at strong 
fields. In the future, we plan to apply this approach for calculations of the transition energies in 
berylliumlike ions that are of current experimental interest [11, 16, 22], 
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